TOPOLOGICAL EXPANSION IN THE CUBIC RANDOM 

MATRIX MODEL 



PAVEL BLEHER AND ALFREDO DEANO 

Abstract. In this paper we study the topological expansion in the 
cubic random matrix model, and we evaluate explicitly the expansion 
coefficients for genus and I. For genus our formula coincides with 
the one of Brezin, Itzykson, Parisi, and Zuber [6]. For the expansion 
coefficients of higher genus we obtain their asymptotic behavior as the 
number of vertices of the associated graphs tends to infinity. Our study 
is based on the Riemann-Hilbert problem, string equations, and the 
Toda equation. 



1. Introduction and statement of the main results 

In this paper we return to the classical work |6j by Brezin, Itzykson, Parisi 
and Zuber, in which, among other things, the authors explicitly calculated 
the coefficients of the topological expansion in the cubic random matrix 
model in genus 0. Our main goal will be to rigorously prove the results of 
[6] and to obtain an explicit formula for the coefficients of the topological 
expansion in genus 1. We will also prove some formulae and asymptotic 
results for the coefficients of the topological expansion in higher genera. 

We consider the random matrix model given by the probability distribu- 
tion 

d m (M) = -l e - w WdM, (1.1) 

on the space of N x N Hermitian matrices M, where 

M 2 

V(M) = — -uM 3 , (1.2) 
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and u > 0. The model is ill-defined because of the divergence at infinity of 
the partition function, 



Z N {u) = J e- N ™^dM. 
The partition function of eigenvalues, 



Z N (u) 



N 



l<j<k<N j=l 



-N[ -j—uz'f 



(1.3) 



dz\ . . . dzN (1.4) 



diverges as well on the real line. To regularize it, we will consider integration 
on a specially chosen contour T in the complex plane: 



Z N (u) 



l<j<k<N 
on which the integral converges. 



N _ N \ ~_±_ uz 3 



dz\ . . . dzjy, (1.5) 




Figure 1. The sectors Sq, Si, S2 and the contours Tq, Ti, T2. 



To choose T, consider the three sectors on the complex plane, 
S = {; 



57T 7tt 1 

z£t: — < arg z < — > 



6 



6 J 



S 2 = [z £ 



7T 



z G C : — < arg z < 
6 



7T 7T 

1 <argz<-- 



}• 



see Fig[TJ Then for any ray 

Re = {z € C : arg z = 6*}, 
lying in the sectors So, Si, and S2, the integral 



z e 



k -N[^-uz" )dz 



Re 



(1.6) 



(1.7) 
(1.8) 
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converges for any k = 0, 1, . . . and any u > 0. Therefore, we will use contours 
consisting of two rays in the sectors So, Si, and #2. More specifically, we 
may consider the three contours on the complex plane, 

ro = Rw U R n /5, Ti = Rn U -R-tt/5, T 2 = R^ 7T / 5 U R n / 5 , (1.9) 

with orientation from (—00) to (00 e 7 ™/ 5 ) on Tq, from (—00) to (00 e -7 ™/ 5 ) 
on Ti, and from (ooe _7r * / ' 5 ) to (ooe ir! '' 5 ) on T2, see Fig. [TJ However, to get 
the topological expansion of the free energy, we will use r = Tq or T = Ti, 
but not r = T2- 

More generally, following Duits and Kuijlaars [13J, it is convenient to 
introduce a linear combination of the contours Tq, Ti. To that end, let us 
fix some a £ C and consider T as a linear combination of Tq, T±, 

T = aT + (1 - a)r x , (1.10) 

in the sense that 

f{z)dz = a [ f(z)dz + {l-a) f f(z)dz, (1.11) 

cf. [1]. With this choice of T = T(a), the integral Z^{u) = Z^{u;a) in 
(|1.5p is convergent for any u > 0. By the Cauchy theorem, we have some 
flexibility in the choice of the contours Tq, T\ within the sectors So, S\, S*2. 

It is easy to see, by differentiation of the integral with respect to u, that 
the partition function Z^{u;a) is analytic for u > 0, and it is infinitely 
differentiable for u > 0. Observe, however, that Z^{u;a) is not analytic at 
u = 0. Indeed, the analytic continuation of Zn(u), u > 0, to u = re~ 3id , 
r > 0, can be obtained by the change of variable Zj = w*e . 



Z N (re- m ; a) = f . . . f \{ ( Wj e ie - w k e i9 f 



l<j<k<N 



N x I .. ,.„.:; 

X 



(1.12) 



Y[e V V d{ Wl e ie )...d{w N e ie ). 

3=1 

By the Cauchy theorem, one can use the same contour T = T(a) for all 9. 
For 9 = ?, the exponential term in the latter formula becomes 

AM -^-+ru>| 
e V 

and the integral in (|1.12p diverges at r = 0. Moreover, 

lim|Zjv(re- 3i *;a)| =00, 9 = ^-, (1.13) 

hence Z^{u;a) is not analytic at u = 0. 
We define the free energy as 

FN ^ a) = W ln Z^a)- (L14) 
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Observe that Zn(0; a) is independent of a. The aim of this paper is to inves- 
tigate the large N asymptotic behavior of the free energy, and, in particular, 
the structure of the different terms that appear in this asymptotic expan- 
sion. We consider u in the interval < u < u c , where u c is the following 
critical value: 

3 l/4 

n c = -j^- (1.15) 
We prove the following results: 

Theorem 1.1. Suppose that a E C is fixed. Then for any < u < u c , the 

free energy F^{u;a) admits an asymptotic expansion in inverse powers of 
N 2 : 

^M~E^. (Lie) 

g=0 

which is uniform in the variable u on any interval [0, u c — e], e > and 
which can be differentiated on [0, u c — e] with respect to u any number of 
times, with a uniform estimate of the error term with respect to u G [0, u c — 
e]. The functions F^ 2a \u) do not depend on a and they admit an analytic 
continuation to the disk \u\ < u c in the complex plane, and if we expand 
them in powers of u, 

00 A2g) 2 j 

**'(«) -£^r-. (»•»> 

then the coefficient f^ 9 ^ is a positive integer number that counts the number 
of 3-valent connected graphs with 2j vertices on a Riemann surface of genus 
9- 

The existence of the -A3 asymptotic expansion for the free energy is known 
in the physical literature since the classical work of Bessis, Itzykson, Zuber 
[7] (see also references therein to the earlier works). A rigorous proof of the 
asymptotic expansion for the free energy for a general polynomial V(M) 
of even degree is given in the paper of Ercolani and McLaughlin [15 1. Our 
proof of Theorem 11.11 follows the work |15j , with some modifications related 
to the fact that the cubic model is defined on the contour T in the complex 
plane and not on the real line. 

It is worth noticing that the coefficients F^ 2 ^ (u) of asymptotic expansion 
(I1.16P in powers of 1/N 2 do not depend on a for all < u < u c . The 
dependence on a arises, however, in the double scaling limit as u — > u c . We 
will return to this question in a subsequent paper. (See also the paper [13] 
of Duits and Kuijlaars, where a similar double scaling limit is studied for 
the quartic random matrix model.) 

In the next theorem we evaluate explicitly the coefficients f$ in the genus 
term F^°\u): 



TOPOLOGICAL EXPANSION IN THE CUBIC RANDOM MATRIX MODEL 5 

Theorem 1.2. The coefficient can be written as 

AO) = 72 J r(f)(2j)! , lgl 

J2] 2 r(i + 3)r(i + i) ! 

and it has the following asymptotic behavior as j — > oo: 

/* m = ^ + °cr'». *, = -£., (1 , 9 , 

where u c is the critical value defined in (|1.15|) . 

From formula ()1.18j) we obtain the first several terms of the Taylor series 
of F(°)(u) at u = 0: 



= 6u 2 + 216u 4 + 13608u 6 + 1119744u 8 + 



540416448 



-u 
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5 (1.20) 
+ 0(u 12 ). 

Our formula (|1.18j) for fjfj coincides with the one obtained in the work 
of Brezin, Itzykson, Parisi and Zuber [6j §4] in the planar diagrams approx- 
imation. 

A more complicated expression, but still explicit, is available when g = 1: 

(2) 

Theorem 1.3. The coefficient can be written explicitly in terms of a 
3F2 hypergeometric function as follows: 

(2) 5-72^T(|)(2 J )! / -j + 1,2,6 . 3\ 

^ ~ 48(3, +2)1^ + 1)1X1 + 1) ^ U "¥ + 1 ' 2 J ' (L21) 

and zi /ias t/ie following asymptotic behavior as j — > 00: 



'« „2ij V' T " U V' "< 48 

Remark 1.4. T/ie 3^2 function in formula (|1 .21 [) can actually be written 
as a linear combination of two 2 i<i functions: 



2,6,-j 3 

3r2 1 k _3j _ 1 ' 9 
' 2 2 z 



2i M -M. _ 1 ' 9 

V 2 3 /' _ (1.23) 

From formula fjl.21f) , the first several terms of the Taylor series ofF( 2 )( It ) 
at u = are 

F^{u) = -u 2 + 189u 4 + 26892u 6 + 4076568u 8 + 3213210384 M io 

2 5 (1.24) 

+ 0(u 12 ). 

In the case of genus g > 1, an explicit formula for the coefficients ffjjj 
becomes complicated, but it is possible to obtain the asymptotic behavior 
of these coefficients when j — > 00: 



6 PAVEL BLEHER AND ALFREDO DEANO 

Theorem 1.5. For any g > 1, the coefficient f^ 9 ^ has the following asymp- 
totic behavior as j — )■ oo: 

fS 9) = K29(2 f — (i + o(r 1/2 )), d-25) 



where 



K 2g = jz A 2g , 1.26 

u 9 



2 

and C 2g satisfies the recurrence relation 

r i / (5 ff -6)(5g-4)cy 2 \ 

- 2 3/2 3 5/4 I 4^ + 54 2^ C 2mC 2m / I 

m+m'=g 
m,m' <g— 1 

(1.27) 

/or <? > 1, mf/i f/ie initial value 

C = _ 2 - 1 /2 3 -7/4 > 

For g = 2 we have that 

8505u 6 „ 8 3905028468u 10 12n 

F (4) = + 2217618u 8 + + 0(u 12 , (1.29) 

2 5 

and the constants C4 and are 

49 • 2 1 / 2 3 3 / 4 7 
= 4ij 6 K = ' . (L30) 

17915904 ' 1440vW V ; 

We note that if we write the recursion (jl.27p as follows: 

C 2g = /z (5g - 6) (5g - 4) C 2g - 2 + v ^ C 2m C 2m ,, (1.31) 

m+m'=g 
m,m'<g— 1 

and if we construct the following generating function (cf. [3]) 

oo 

y(t) = Y t C 2g t^ L , (1.32) 

9=0 

then y(t) satisfies the Painleve I differential equation 

/'(O^-g, (1.33, 

We can make the change of variables t = —ct and n = Ay to bring it to the 
standard form 

u"( T ) = 6u 2 {t) + t, (1.34) 

see p]. Explicitly, 

c = 2 - 3 / 5 , A = 2 3 / 10 3 5 / 4 . (1.35) 

The proof of Theorem 11.11 relies on the fact that the equilibrium measure 
for this problem is supported by a single interval [a, b] on the real axis, 
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provided < u < u c . This will be shown in Section [2j As a consequence, we 
will be able to apply the technique developed in [15] to the proof of Theorem 

o 

The precise structure of the coefficients F^ 2 ^ (u) of the asymptotic expan- 
sion of the free energy will be analyzed using the following steps: 

• First we apply the Deift-Zhou nonlinear steepest descent method to 
the associated Riemann-Hilbert problem, in order to obtain the large 
TV" asymptotic expansion of the recurrence coefficients 7^ and /3 n of 
the corresponding orthogonal polynomials P n (z), when the index n is 
of the order of N. The coefficients in this asymptotic expansion will 
be found by using the string equations, which are nonlinear algebraic 
relations for 7^ and /3 n . Simultaneously, this proves the existence of 
the orthogonal polynomials P n (z) when the index n is of the order 
of ./V and iV is large enough. 

• Then we make a change of variable to reduce the cubic polynomial 
V(M) = ^ - uMS to V(M) = tM - and we derive simple 
formulae for the partition function and the recurrence coefficients 
under this change of variable. The change of variable allows us to 
use the Toda equation that connects the second derivative of the 
free energy Fn with respect to the parameter t to the recurrence 
coefficient 7^. 

• By integrating the large N asymptotic expansion for 7^ term by 
term, we obtain detailed information about the large N asymptotic 
expansion of the free energy Fn (u) , and in particular the form of the 
different F^ 29 \u) terms. Note that in this derivation we only need 
the recurrence coefficients 7^ and /3jv for large values of N, whose ex- 
istence and asymptotic behavior are proved via the Riemann-Hilbert 
analysis before. 

We note that an alternative method to derive formulae for F^ 2 ^ was pro- 
posed in [2], using the so called loop equations (see [16] for a rigorous deriva- 
tion of the loop equations). In [2], explicit expressions for F® and F^ arc 
given, in terms of elementary functions that depend on the endpoints of the 
support of the equilibrium measure associated with the corresponding po- 
tential. In [HI [19], a general heuristic formula for computing F^ 9 ' is presented 
for arbitrary potentials and multicut cases. 

We would like also to bring attention to the very recent paper [18] of Er- 
colani and Pierce, which appeared after the present paper had been posted 
on arXiv, and which cited the present paper. In [18] a different approach, 
based on the so called difference string equations, is developed for the calcu- 
lation of the terms of the topological expansion in the random matrix model 
with a general potential V(M). 

The rest of the paper is organized as follows: in Sections [2] and [3] we 
analyze the equilibrium measure and the corresponding Riemann-Hilbert 
problem for orthogonal polynomials with respect to the cubic- type weight. 
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In Section H] we apply the string equations to the large N expansion of the 
recurrence coefficients 7^ and f3 n . This is used in Section [5] to obtain the 
large N expansion of the free energy Fn(u), and this leads to the proof of 
Theorems 11.21 and Theorem 1 1.3 1 in Section [6j In Section [7] we prove Theorem 
11.51 expanding the recurrence coefficients around the critical point u = u c 
and selecting the terms that give the dominant behavior in the free energy. 
Finally, Section is devoted to the interpretation of the topological expan- 
sion in terms of connected graphs embedded in closed Riemann surfaces. 



2. Equilibrium measure and the Riemann-Hilbert analysis 

2.1. Construction. Let us denote by g(s) the density of the equilibrium 
measure for this problem, that we assume supported on a certain curve 
J C C with endpoints z = a and z = b. We consider J oriented from a to 
b, and we take the +- side on the left of J and the — side on the right of J, 
following the standard convention. 
The resolvent 

= f g(s)ds 
J J z-s 

is analytic off J, and it satisfies the Euler-Lagrange equation on J, 

u+(z)+U-(z) = V(z), zeJ, (2.2) 
and the asymptotics at infinity, 

u(z) = - + 0(z~ 2 ), z-^oc. (2.3) 

To solve (12.21). we write 



u J (z) = -V'(z)--^R(7)h(z), (2.4) 

where h(z) is an analytic function and 

R(z) = (z- a)(z-b), (2.5) 

since we are assuming that we are in the one cut case. We take the principal 
sheet for \J R(z), with a cut on J. Due to the Plemelj formula, we have 

OJ-(z) - oj + (z) = 27rig(z) = J R(z)h{z) g{z) = — J R(z)h(z). 

We easily deduce from fj2.4j) that h(z) is a polynomial of degree 1. Moreover, 
if we write 

a = x — y, b = x + y (2-6) 
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and h(z) = A(z — zq), then identifying the asymptotics at infinity on both 
sides of (|2.4j) . we find 

A = -3u, 
1 

zo = o X ' 

Su 

2 2 2x ( 2 - 7 ) 
2x 2 + y 2 = — , 
6u 

y 2 

— (1 — 6ux) = 1. 



Thus, 



co(z) = —y/(z — a){z — b) (1 — 3nz — 3ux) (2. 



and therefore 



q(z) = ^-V( z -o){z- b) (1 - 3uz - 3ux). (2.9) 
27^^ 



It follows that g(z) 2 has a double zero at 

1 

z 0= o X 

6u 

and two simple roots located at z\ } 2 = a, b. 
Next we will prove the following facts: 

• For < u < u c , where u c is given by (|1.15p . the equilibrium measure 
is supported on an interval J = [a, b] of the real axis. The endpoints 
of this interval are analytic functions of the parameter u. 

• For < u < u c , we can extend the interval J to an unbounded 
contour Tq = 71 U J U 72 in the upper half-plane, 

C+ = {z G C : Imz > 0}, 

in such a way that 

M z )>® 2G71, 
(j> 2 (z) > z £ 72, 



(2.10) 



where 



1 



i(*) = o / VR(s)h(s)ds 



(2-11) 



02(2) = - / v ^R(s)h(s)ds. 



Additionally, as |#| — ^ 00 on 71 and 72, it is possible to estimate the 
growth of Re <f>i (z) and Ke<p2(z): 

Re^(z) = -^ + 0(z 2 ), z = 1,2. (2.12) 
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This last result is relevant in the steepest descent method applied 
to the Riemann-Hilbert problem, in order to show that the jump 
matrices outside of J tend to the identity when N — > oo. 

2.2. Support of the equilibrium measure. Combining the last two equa- 
tions in ()2.7p . we obtain the following cubic equation for the variable 
function of the parameter u: 

I8u 2 x 3 - 9ux 2 + x - 6u = 0. (2.13) 

If u > is small, then this equation has the three solutions, 

x 1 = 6n + 0(n 2 ), x 2 = ^- + 0(l), x 3 = ^- + 0(l). (2.14) 

bu 6u 

In what follows we will be interested in the first root, x = x%, which is the 

solution that remains bounded when u is small. More terms in the expansion 

of x are 

x = xi{u) = Qu + 324-u 3 + 31104n 5 + 0(u 7 ). (2.15) 
From the last equation in (|2,7p we find 

y = 2 + 36u 2 + 2916u 4 + 0(u 6 ). (2.16) 

The support of the equilibrium measure is the interval [a, b] = [x — y,x + y]. 
As u — > 0, it converges to [—2,2], which is the support of the equilibrium 
measure for GUE. 

The series that we obtain in (|2.15j) can be proved to be convergent for 
small values of u: 

Proposition 2.1. The series (12.150 is convergent for \u\ < u c , where u c is 
given by (I1.15p . 

Proof. Since by (|2.14p the root x = xi(u) is isolated for small u, xi(u) 
is analytic at u = 0. The discriminant of cubic equation (|2.13p is A = 
9u 2 (l — 34992m ), and it vanishes at the critical value u = u c . Hence it does 
not vanish for < |u| < u c , and the root x\{u) remains isolated. This proves 
that x\(u) is analytic in the disk \u\ < u c , and hence the Taylor series for 
x\{u) is convergent in this disk. □ 

It follows from the last equation of (|2.7p that a similar result is valid for 
the variable y = y{u) > 0. Namely, y(u) is obviously analytic at u = 0. 
Suppose that 1 — 6ux\(u) = 0, then from equation (|2,13p we obtain that 
u = 0, which contradicts 1 — Qux\(u) = 0. Therefore, 

1 — Qux\(u) ^ 0, \u\ < u c , (2-17) 

hence y{u) is analytic in the disk \u\ < u c . This shows that for real < u < 
u c , both x = x\(u) and y(u) > are real, so the support of the equilibrium 
measure is the interval [a, b] on the real axis. Moreover, the endpoints z = a 
and z = b are analytic functions of the parameter u. 

In order to prove that we are in the one-cut regular case, we also need to 
show that the double root zq lies outside [a, b]: 
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Proposition 2.2. For < u < u c , we have that zq > b, so that the double 
root is outside the interval [a, b] . 

Proof. By the second equation in (|2.7p . 

1 

zq = - x ->■ +oo, 

6u 

as u — > + , so zq > b for small u > 0. Suppose that for some < u < u c we 
have zo = 6. Then 

hence 

2x) = y 2 = 

J 1 — 6ux 

and 

(1 -6uxf = 36u 2 . 

Denote v = ux. Then the latter equation and (|2.13j) give two cubic equations 
on v. 

(1 - 6vf - 36u 2 = 0, 18v 3 -9v 2 + v- 6n 2 = 0. 
By applying the Euclidean algorithm to these two cubic equations, we obtain 
that it 2 (l — 34992n 4 ) = 0, which is not true for < u < u c . Therefore, zq > b 
for < u < u c . □ 

When u = u c , the support of the equilibrium measure will be the interval 

[a, b] = [3 3 / 4 - 3 5 / 4 , 3 3 / 4 + 3 1 / 4 ], (2.18) 

and the double root of the function g(z) 2 is located at 

Zc = 3 3 / 4 + 3 1 / 4 , (2.19) 

which coincides with the right endpoint of the support of the equilibrium 
measure. 

We prove now that there exists an extension of the interval J to a curve 
r in C with right properties for the Riemann-Hilbert analysis. 

Proposition 2.3. There exist two curves 71 and 72 in C such that condi- 
tions (|2.10p hold true. 

Proof. If b < z < zq then (z — a)(z — b) > and z — zq < again, so 
Re^2(z) > and we can take the segment (b,Zo) as part of 72. However, 
we cannot take the whole interval (6,00) as 72, since Re <fo(z) < if z > zq. 
Nevertheless, if we define the quadratic differential —Q{z)dz 2 , where 

Q{z) = \R(z)h 2 (z), (2.20) 

then it follows from the general theory (see \25\ §8.2] or [271 Theorem 7.1]) 
that, since z = zq is a double zero of —Q(z)dz 2 , then there are exactly four 
trajectories that emanate from z = Zq, along which ImQ(z) = 0. Locally, 
these trajectories leave the point z = zq with angles ki:/2, < k < 3. 
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■ 










a 


b 


^0 



Figure 2. The contour f = 71 U [a, b] U 72. 



We choose the trajectory that leaves z = Zq with angle n/2 into the upper 
half plane. This trajectory must go to 00, since ^2(^0) > and faiz) is 
increasing on 72, so it cannot go back to the real axis. Now, 



< (f> 2 (z) - Mzo) = \ I VW)Ks)d& 

^ J ZD 

u(s) \ ds 



zo 

z 'V'(s) 



Jz \ 



(2.21) 



V(z) - V(z ) 



uz 



+ 0{\n\z\) = -^ + 0(z 2 ) 



when I z \ — > 00. Since (j>2 (z) is real for z G 72, we have Re<fo(z) = <j>2 (z), and 

^3 

(2.22) 



uz 



ReMz) = -^ + 0{z 2 ) > 0. 



This also shows that 72 goes to infinity with angle 7r/3. 

We can take 71 = (—00, a), since for z < a we have (z — a)(z — b) > and 
z — zo < 0, so 

IVWMz) = ^y/( z - a )(z-b)(z - z ) < 0, (2.23) 

and hence Re <j>i(z) > for z < a, after integrating. Note that on (—00, a) 
we need to take the square root with negative sign. On (—00, a), with a 
similar computation as before, we have 

.3 

(2.24) 



uz 



ReMz)=Mz) = -^-+0(z 2 ), 



and the result follows. 



□ 
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2.3. The Riemann Hilbert analysis. We will consider monic orthogonal 
polynomials (OPs) P n {z) = z n + . . . with respect to the weight function 

w(z) = e- NV{z \ (2.25) 

so that 

J P n (z)z k w(z)dz = 0, fe = 0,l,...,n-l. (2.26) 

The existence of the OPs is not obvious, since the contour T is complex, but 
the Riemann-Hilbert (RH) analysis provides us, among other properties, 
with the existence and uniqueness of P n {z) for large enough values of n. 
The RH analysis also gives us various properties of the OPs. In particular, 
we will see that the OPs satisfy the three-term recurrence relation, 

zP n (z) = P n+l {z) + p n P n {z)+ 1 2 n P n . l {z), (2.27) 

and that the recurrence coefficients fi n , 7^ satisfy string and deformation 
equations. 

It is known from the work of Fokas, Its and Kitaev [20], that the family of 
orthogonal polynomials can be characterized as the solution of the following 
2x2 Riemann-Hilbert problem (RHP): 

Find a 2 x 2 matrix-valued Y n : C \ (r U Fx) -> C 2x2 such that 

• Y n {z) is analytic for z G C \ (ro U T±) and the limits 

lim Y n (s) = Y n± {z) 

S->2±0 

exist as s approaches z from the ±-side of To U Ti. As usual, we 
assume that the +- side is on the left of an oriented contour. 

• For z G ToUTi, 

Y n+ {z) = Y n _{z) (J a (^~ NV{z) ^ f (2.28) 

where 

a(z) = a X r Q {z) + (l-a)xr 1 (z), (2.29) 
and xa{z) denotes the characteristic function of the set A. 



As z — > 00, 



V '" (:) ~l / + C, ^))(^ A' l2 ' :iU) 



We will call n the degree of the RHP. This RHP has a unique solution if and 
only if the monic polynomial P n (z), orthogonal with respect to the weight 
function w(z)e~ NV ( z \ uniquely exists (see the next section). If additionally 
P n _i(z) uniquely exists, then the solution of the Riemann-Hilbert problem 
is given by: 

VM _( (CP nW )(z) \ 

n[Z) " \-&Pn-l(z) -j^(CP n ^ W )(z) " (2 - 31) 
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where 




ds 



(2.32) 



is the Cauchy transform on T, and the coefficient h n —\ is defined as 



The nonlinear steepest descent method, due to Deift and Zhou (see 
and also [12]), consists of a series of explicit and invertible transformations 
of this original Riemann-Hilbert problem, with the objective of arriving at 
a transformed problem for a 2 x 2 matrix R(z), where the matrix jumps 
tend to the identity uniformly in C. These steps are briefly as follows: 

• Y n i— >■ T n , normalization at oo, using the equilibrium measure on J. 

• T n I—?- S n , opening of lenses around J to convert the highly oscillatory 
jump on J into a combination of exponentially close to identity and 
constant matrix jumps. 

• Sn i — ^ R n , construction of global and local parametrices, which are 
(respectively) approximations near the interval J and close to the 
endpoints of it. In this case, the local approximation near the end- 
points z = a and z = b are given in terms of Airy functions. 

We refer the reader to the works [12], [10], [15] for the details. In this 
case, the Riemann-Hilbert analysis gives the following result: 

Theorem 2.1. For any 5 > 0, there exist Nq(5) > and e(5) > such that 
for any u G [0,u c — 8] and any N > Nq(5), a solution Y n (z) to the RHP 
(|2.28p -( i2.30p exists for n in the interval 



One of the key results also given by the Riemann-Hilbert analysis is that 
the recurrence coefficients 7^ and j3 n admit an asymptotic expansion in 
inverse powers of N for large N. We will return to these expansions later, in 
Section HI but first we need to discuss various properties of the orthogonal 
polynomials, which follow from the RHP. 



The orthogonality conditions (|2.26[) are equivalent to a linear system of 
equations for the coefficients of the orthogonal polynomial 




(2.33) 




(2.34) 



3. Orthogonal polynomials P n (z) 



P n (z) 



z n + a n -iz 



n— 1 



+ . . . + a\z + a . 



(3.1) 



Namely, if 




(3.2) 



then 



n-l 




k = 0, 1, . . . , n — 1. 



(3.3) 



3=0 



TOPOLOGICAL EXPANSION IN THE CUBIC RANDOM MATRIX MODEL 15 



Therefore, the orthogonal polynomial P n (z) exists and it is unique if and 
only if linear system (13.3H is nondegenerate, so that 



In this case, 



where 



D n (z) 



C Cl 
Cl c 2 

Cn— 1 C n 
Pn(z) = 



C Ci 
Cl C 2 



Cn— 1 c n 
1 2 



Cn-1 
C2n-2 



^0. 



(3.4) 



n-l 



(3.5) 



Cn— 1 C n 
C2n-2 C2n-1 



2 



n-l 



Z 



(3.6) 



see [28]. In addition, we have 

Proposition 3.1. If there exists a unique orthogonal polynomial P n (z) then 

D n = h n D n - 1 , (3.7) 

where 



h r , 



P n (z) 2 w(z)dz. 



(3.8) 



For a proof of this and subsequent propositions of this section see Appen- 
dix [A] at the end of the paper. Another useful formula for D n —\ is 



A.-1 = \ I • • • / A(z) 2 f[ e~ NV ^d Zl . . . dz n , 
n - J r Jt j= i 



(3.9) 
(3.10) 



where 

n fe--^) 2 - 

see [28]. 

Now we will relate the existence of a solution to RHP (|2.28|) - (|2.30p to the 
existence and the uniqueness of orthogonal polynomials. 

Proposition 3.2. Suppose that RHP ([23g]) - (pOD|) has a solution Y n (z). 
Then there exists a unique orthogonal polynomial P n {z). 

The following proposition proves the existence of the three term recur- 
rence relation: 

Proposition 3.3. Suppose that RHP (j2728|) - (f230|) has solutions Y n -x{z), 
Y n (z), and Y n+ i(z) for the degrees n — 1, n, andn + 1, respectively. Then the 
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orthogonal polynomials P n -\{z), P n {z), and P n+ \{z), which uniquely exist 
by Proposition 1 3. S\ satisfy the three term recurrence relation, 

zP n (z) = P n+1 (z) + PnP n (z) + l 2 n Pn-l{z). (3.11) 

Observe that under the conditions of Proposition 13,31 D n / 0, D n -\ 7^ 0, 
and Dn-2 7^ and hence, by (|3.7p . 

K = -^± 0, V-i = #^ / 0. (3.12) 

From (I3.1ip we obtain, by multiplying by P n —\{z) and integrating with 
respect to w(z), that 

7n = t^-^O. (3.13) 

Note that the proof of Propositions I3.HI3.3I does not use specific form 
(|2.25p of the weight w(z), and therefore they hold in a much more general 
situation. In the two subsequent propositions we assume that the weight 
w(z) has form (|2.25p . 

Proposition 3.4. Suppose that RHP (pT28|) - (f2T30l) has solutions Yj(z) for 
the degrees j = n, n ± 1, n ± 2. T/ien £/ie string equations hold: 

3«(7n+l + Pn + 7n) = An 

n (3.14) 
7 2(l-3n(/3 n + /3„_ 1 )) = -. 

By Theorem 12, H for any 5 > 0, Propositions I3.1H3.4I are applied for u G 
[0, u c — 8], provided {n, n±l,n±2} C In = In (8). Observe that In (8) does 
not depend on u. 

Proposition 3.5. For any 5 > 0, if {n,n ± l,n ± 2} C In(S), then the 
coefficients (3 n and 7^ are C°° -functions of u for u G [0, u c — 5] and i/iej/ are 
analytic functions of u for u £ (0, u c — 5]. 

Remark. In what follows we use the recurrence coefficient 7^, which is 
uniquely defined, and we will not use 7n , which is defined up to a sign. If 
needed, the sign of 7„ can be uniquely defined by the analytic continuation 
of 7 n in u starting from the value 7n (0) = yn/N. Note, that, in general, 
the recurrence coefficients j3 n and 7^ are complex valued. 

4. String equations and the large N expansion of the 

RECURRENCE COEFFICIENTS 

4.1. The large N expansion of the recurrence coefficients. From the 
Riemann-Hilbert analysis, we know that 7^ and (3 n can be expanded in 
inverse powers of N, as functions of the parameter s = n/N. Moreover, in 
the expansion of 7^ the odd coefficients vanish, and we have an expansion in 
inverse powers of N~ 2 . Similarly, for /3 n we also have an expansion in inverse 
powers of N~ 2 , if we consider (3 n as a function of the shifted parameter 
s = -ft + 277. Namely, we have the following result: 
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Theorem 4.1. (See [5]). For any 5 > 0, ifn C In(5/2), then the coefficients 
Pn, 7^ can be expanded in uniform asymptotic series, 



fc=0 
oo 

E 

k=0 





( n 




ViV 




( n 







(4.1) 



where the functions g2k{s, u), b2k{s, u), k = 0, 1, . . do not depend on n and 
N , and they are C 00 -smooth in s and u on the set {0 < u < u c — 5, \s — 1| < 
e(5)/2}. 

The proof of this theorem in [5] is based on the string equations and 
it is applied to the current case without any changes. To find the func- 
tions g2k(s, u), &2fc(sj u) iteratively, we substitute expansions (|4.ip into string 
equations (I3.14j) and equate coefficients at different powers of N~ 2 . Let us 
consider the first equation in (|3.14p . If we take s = jj + 
first equation in (|4.ip we have the asymptotic expansions 



+ 77X7, then from the 



7n ~ 



:92k s 



1 

2iV' 



j\r 2fc 

fc=0 

°° 1 / 1 \ 

n+^J2j^k92k[s + —,u) , 
u—n \ / 



n 1 

N + 2N 



(4.2) 



We substitute these expansions into the first equation in (|3.14p and expand 
functions the g^k into Taylor series at s: 



3u 



2 EE 

k=Q j=0 



92k\s,u) 
{2j)\2 2 i N 2k + 2 i 



+ E 



\k=0 



hk(s,u) 

N 2k 



hk(s,u) ri 



k=0 



where 



(2i)/ \ 
92k (s,u) 



d 2] 92k{s-,u) 
ds 2 i 



(4.4) 



By equating the coefficients at the same powers of N~ 2 on both sides, we 
obtain a series of equations for the functions g2k, 02k- 

Similarly, if we take s = -K>, then from the second equation in (|4,ip we 
obtain that 



E 

\k=0 



92k(S,U) 
N 2k 



i-HEE 



fc=0 j=0 



(2j)!2 2 J'7V 2fc + 2 i I 



(4.5) 



Let us analyze the equations obtained when we equate the coefficients at 
the same powers of N~ 2 on both sides of equations (|4.3p . (|4.5p . 
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4.2. Zeroth order terms in the recurrence coefficients. At the zeroth 
order of N~ 2 we obtain the system of equations, 

3u[2g (s, u) + b (s, u) 2 } = b (s, u) , (4.6) 
g (s,u)[l - 6ub (s,u)] = s. 
These equations can be solved to yield 

72u 2 g (s, uf - go(s, u) 2 + s 2 = 0, 

, / \ 9o(s,u) - s (4.7) 

O0{S,U) = — -, r- , 

6ug (s,u) 

or alternatively, 

186o(s, u) 3 u 2 — 9bo(s, u) 2 u + 6q(s, u) — 6su = 0, 

8 (4.8) 

^ ^ S ' U > = 1 eTlTl V • 

1 — oubo(s, u) 

Consider s in a small neighborhood of the point s = 1 on the complex plane, 
say, \s — 1| < c <C 1. Then the cubic equation (|4.7p admits three different 
solutions, which behave as follows when u — > 0: 

g i i\s,u) = s + 0{u 2 ) 

g^(s,u) = -s + 0(u 2 ) (4.9) 
g^(s,u) = -^ + 0(l). 



The solution we are interested in is go(s,u) = g$ (s,u). Since g^\s,u) is 
an isolated solution of the cubic equation, it is analytic at u = 0. In fact it 
is analytic in the disk \u\ < u c s~ l l 2 on the complex plane. In the following 
proposition we explicitly evaluate the Taylor series of go(s,u) = g^\s,u) at 
u = 0. 

Proposition 4.1. The Taylor series 

oo 

g Q (s,u) = Y,a 2j (s)u 2 i ( 4 .1 ) 
i=o 

converges in the disk \u\ < UcS^ 1 ! 2 and its coefficients are 

72 i c i+1 T( 3j+1 ) 
a 2 i(s) = 2 ' , j > 0. (4.11) 

3K ' 2r(j + i)r(^)' y ' 

Proof. We make the change of variables, 

w = su 2 , v = u 2 go(s,u). (4-12) 
Then v = v(w) satisfies the cubic equation 

72v 3 - v 2 + w 2 = 0, (4.13) 
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and we are looking for the root v = v(w) such that v(0) = and v(w) > 
for small w > 0. Observe that v(w) admits an expansion in powers of w, 

oo 

v{w) = y]cjW*. (4-14) 

j=i 

Using Cauchy integral formula, if 7 denotes a smooth closed contour around 
w = 0, we have 

Cj = — l V -^dw. (4.15) 

Since v(Q) = 0, we make a change of variables 

1 / v dw 
Cj ~ 2iri w(vy +1 dv 

We have 



dv. (4.16) 



dw 1 - W8v 



w(v) = -VI - 72v, — = , (4.17) 
W s dv VI - 72v V ; 



so 



Cj 2m 



v~ j (l - 108u)(l - v)~~ l dv. (4.18) 
We expand the binomial series: 

00 /— i- — 1\ 

(1 - 108v)(l - v)-^ 1 = (l-108w) V 2 , (-l) fc (72?;)* 

fc=cA fc ' 

= (l-108,)f;( fc + ")(72^, 

and we pick up the residues, that correspond to k = j — 1 in the first term 
and k = j — 2 in the second one. Then 



(4.19) 



Cj = 72^ 
Now, 



f-l\ 3/f-2 



j-iy 2 Vi-2 



2r(i)r(| + i) 



(4.20) 



00 00 

Sr. S 



9o(s, u) = —v(w) = c o w3 = s ^2 c j+i wJ 

W 3=1 3=0 



72V +1 r( 



(4.21) 



2 



2r(j + i)r(i±5) 



which gives the expression for the coefficients a%j(s). □ 
It is easy to generate several coefficients using general formula (|4.1ip : 
g (s, u) =s + 36sV + 3240sV + 373248s 4 u 6 + 48498912s 5 u 8 
+ 0[u 1 . 

As a consequence of Proposition 14.11 we have 



( 4 - 22 ) 
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Corollary 4.2. The coefficient bo(s, u) is an analytic function of u for \u\ < 
u c s- l /\ " 

Proof. From the second equation in (|4.7|) and formula (|4.22|) it is clear that 
bo(s,u) is analytic at u = 0. In addition, go(s,u) / for |u| < u c s -1 / 2 , 
because if go(s,u) = then by the first equation in (14. 7h s = 0, but we 
assume that |s — 1| < c< 1. Therefore, the analyticity of go(s,u) implies 
the analyticity of bo(s,u) in the disk |u| < u c s~ 1 / 2 . □ 

To obtain the asymptotics of 7^ and (3n as N — > 00, we set n = N in 
(|4TTjh so that 



7 ^~<7o(l,n) + ^^ 52fc (l,n), (4.23) 
k=l 

P» ~ b ° (* + ^ u ) + E j^S 6 ^ (l + u) . (4.24) 
In the sequel, we denote 

02fc(«) =#2fc(l ! «)> ^2fe(^) = &2fc(l,«), 

for brevity. When s = 1, the cubic equation (|4.7p for go(u) becomes 

72u 2 g (u) 3 - g (u) 2 + 1 = 0, (4.25) 
and the solution admits the power series expansion 

valid when |u| < u c . Formula (|4.22p reduces for s = 1 to 

5o (u) = 1 + 36u 2 + 3240u 4 + 373248m 6 + 48498912tt 8 + 0(u 10 ). (4.27) 

4.3. Analysis of the higher order terms. An important consequence of 
string equations (|4.3|) and (|4.5j) is that higher order terms g2k = 92k{s,u) 
and &2fc = b2k(s,u), k > 1, can be expressed in terms of <7o(sj^)i &o(s,u), 
and their derivatives with respect to s. More precisely: 

• We can solve for g^h and &2fc i n terms of the previous coefficients 
through a linear system of two equations. 

• The analysis of the determinant of this linear system shows that 
higher order terms have the same singularity u = ±u c as the ini- 
tial terms go and bo, and that no other singularities appear when 
increasing k. In this respect, this property is analogous to the one 
presented in p3], [15] and [T7] for general even potentials V(M). 

• Although an explicit formula for (^fc and ^2fc is rather complex when 
k > 2, in order to determine the large j asymptotic behavior of the 
coefficients /o? in the expansion of the free energy, it is enough to 
evaluate the leading order term at the singular points u = ±u c , and 
this is possible to extract from the string equations. 
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For convenience, we make the change of variable 



z = £ (4.28) 
u 



so that the potential V(z) becomes 



^) = i(y-^ 3 )- ( 13)) 



With this change of variable, we introduce a new family of orthogonal poly- 
nomials 

P n (0 = u n P n (£) , (4.30) 



and the corresponding recurrence coefficients are 

/3 n = uf3 n , = u^l (4.31) 

Now the string equations (|3.14j) read 

my 2 ( 4 - 32 ) 

t(i -m +/?„-!)) = — . 

We will use the same scaled parameter as before, 

w = su 2 . (4.33) 
From (I4.ip and (j4.32f) we have an asymptotic expansion for these new coef- 

2 

flcients in powers of jr: 

00 ,4k /^„,2- 



2 ST~^ u I nu 

7n ~ 2^ 92k 



N 2k y I 7 ' 

fc=0 

00 4A- / 2 2 \ 

U f / TO U \ 



fc=0 

where 



(4.34) 



92k(w) = u- 4k+2 g 2k (s, u), b 2k (w) = u-* k+1 b 2k (s, u). (4.35) 
For s = 1 this reduces to 

52fe(n 2 ) = u- 4k+2 g 2k (u), b 2k (u 2 ) = u- ik+1 b 2k (u), (4.36) 
and in particular, when k = 0, we have 

g (w) = u 2 g (u), b (w) = ub (u). (4.37) 
It follows from (|4TTj) that 

SoW = £^L^, (4. 38 ) 

and also that go(w) satisfies the following cubic equation: 

72g 3 (w)-g 2 (w)+w 2 = 0, (4.39) 
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which is identical to ()4.13p . The critical value now becomes 

Wc = ul = (4.40) 
c 324 v ) 

By (|4.14|) and (|4.20p . the function go(w) is analytic at w = 0, and it has the 
following Taylor expansion: 

g ( W )=\^^l L . 4.41 

^ 2r(i)r(| + i) 

From equations fj4.3j) and (|4.5p we obtain the string equations for g 2 k{w) 
and b 2 k(w): 

f,^ (.„) 

6 E 7|W + 3 E ^W^'W=^H (4-42) 

m+j=k m+m'=fc 

and 

52fcH-6 2^ (2jjtift ' " ( ^ 

The advantage of equations (|4.38j) — (|4.43|) is that they do not contain the 
parameter u anymore: it is hidden in the argument w = su 2 . 

By solving these equations for g 2 k(w) an d b2k(w), we obtain the following 
linear system of equations for k > 1: 

6g2k(w) + (660 (w) - l)b 2k {w) 

m+j=k m+m'=k 
m<k—l m,m'<k—l 

and 

(1 - 6&oM)02k(«O - Qgo{w)b 2k (w) =6 ^ (2?)!2 2 i 

m+m'+j=k 
m,m' <k— 1 

(4.45) 

The determinant of this system is 

D(w) = -36g (w) + (1 - 6&oO)) 2 - (4.46) 
By using equations f|4.38j) . f|4.39j) . it can be simplified to 

D(w) = 1 - 108<7o(w). (4.47) 
This leads to the following result: 

Proposition 4.3. For k > 1 f/ie functions g 2 k{w) an d b 2 k(w) are analytic 
functions of w in the disk \ w\ < w c in the complex plane. 
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Proof. The proof of this proposition is by induction in k. To this end, we 
need to show that the determinant D{w) does not vanish for \w\ < w c . 

Suppose D(w) = for some w such that \w\ < w c . Then by (14.471) . we 
have 

9oH = ^- (4-48) 
By substituting this into the cubic equation (J4.39P we obtain that 

hence D{w) / inside the disk \w\ < w c . 

The proposition now follows, since it is clear that the only singularities 
of the higher order terms for \w\ < w c are those inherited from go(w) and 
bo(w), but go(w) and bo(w) are analytic in the disk \w\ < w c . □ 

As a consequence of Proposition 14.31 we have the following result: 

Corollary 4.4. The functions g 2k (s,u) and b 2k (s, u ) are analytic in the disk 
\u\ < u c s~ l l 2 in the complex plane, and g2k(s,u) is even in u and b2k(s,u) 
is odd in u. 

Proof. Analyticity follows from Proposition 14.31 and the change of variable 
w = su 2 , and the parity, from rewriting equation (|4.35[) : 

g 2k (s, u) = u Ak - 2 g 2k {w), b 2k (s, u) = u ik ~ l b 2k (w). (4.50) 

□ 

5. The large N expansion of the free energy Fn 

5.1. Another change of variable and an analysis of the auxiliary 
potential. To evaluate the large N expansion of the free energy Fn, we 
will use the Toda equation, which relates Fn to the recurrence coefficient 
7^. To this end we will consider a different form of the partition function. 
Under the change of variable 

2f = (3t*)-V3 C + * (5.1) 
bu 

where we assume that u > and (3-u) -1 / 3 > 0, we have, by straightforward 
algebra, that 

17 " ^ 3 " 77^-2 ="V+*C, (5-2) 



where 



Let us denote 



2 " 108u 2 3 
1 



(5.3) 



4(3n) 4 / 3 ' 

r N 

2 "= - II & - C*) 2 II e~ NViCj) d(i ■ ■ ■ d( N , (5.4) 
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where 

c 3 



V(()=-f + t(, (5.5) 



and r is the image of T under the change of variable (|5.ip . Note that 
this cubic model was used in [9] with t = 0, in the context of the complex 
Gaussian quadrature of integrals with high order stationary points. The 
corresponding free energy is 

F N = j^\nZ N . (5.6) 

Formula (|5,2p implies that 

F N = ^ + ^ + F N = 2 -f^-hm + FN . (5 . 7) 
108u 2 3 3 4 v ' 

By Theorem II. 1\ for any e > 0, F^(u) admits a uniform asymptotic ex- 
pansion for < u < u c — e, hence for any e > 0, F^(t) admits a uniform 
asymptotic expansion as well, for t c + e < t < oo: 

k=0 



where 



^ = lob + ^ + F<0) (») = - ^ + F,0, <")- < 5 - 9 > 

and for fc > 1, 

pm {t) = F m {u)j u = — i^. (5.io) 

From (I1.15P and (15. 3p we obtain that the new critical value is 

t c = , — ^ = 3- 2~ 2/3 . (5.11) 

4(3u c ) 4 / 3 V ' 



Note that with the change of variable (|5.ip . the family of orthogonal 
polynomials is scaled and shifted. If we write z = cC, + d, then: 

P n {z) = P n (c( + d) = c n P n (C), (5.12) 

so that P n (C) is again monic. The coefficients of the recurrence relation are 
modified in the following way: 

a Pn — d 2 In i q\ 

Pn = , In = ~2- (5.16) 



In our case, according to (|4.4U|) . c = (3u) 1//3 and d = J7, so 



7* = (3^) 2/3 7^ 



7-, 



2 



2VT 

1 \ I / .. 1 

6u 



K = (*,)>/» ( A - i-) = ^ (a. 



(5.14) 
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An important fact in our analysis will be that with respect to the variable 
t we have the following Toda equation (see, e.g., [5]): 

d 2 F N _ 2 

(5-15) 

The usual derivation of the Toda equation assumes the existence of the 
orthogonal polynomials P n (C) f° r n = 0,1, . . . ,N. In our case we know the 
existence of the orthogonal polynomials P n {z), and hence P n (C) only for 
n G In- We will prove that, nevertheless, the Toda equation is valid in our 
case. Namely, we will prove the following proposition: 

Proposition 5.1. Toda equation (|5. 15[) holds for any t > t c . 

For a proof of this proposition see Appendix [B] in the end of the paper. 
The idea now is to integrate twice equation (I5.15D in the variable t, in 
order to obtain Ftf, and then compute F/v using equation (|5.7p . From the 
asymptotic expansion (|4.23p . we have 

2 92k(t) 1 92k{ u ) 1 

2.^ = ^/1 2. -jvsr. u = mW l - (5 - 16) 

k=0 k=0 v ' 

Integrating this expression twice should give an expansion for Fn, and then 
we can compute Fn using (|5.7p . However, the problem now is that we 
have to justify that the term-by-term integration of the large N asymptotic 
expansion of over an unbounded interval in the variable t is permissible. 
This is the content of the following proposition: 

Proposition 5.1. We have that 

^)^-^ + /:/;(^)-^-^)^ (-> 

and for any k > 1, 

p(2k ){t)= j r ~ g2k{a)d(TdT , (5 . 18 ) 

J 00 J 00 

Proof. Let us prove first that for any k > 0, 

d 2 F^ 2k \t) . , , , 
-^- L =92k{t). (5.19) 

Consider the difference operator 

&hf(t) = f{t + h)- 2f{t) + f(t - h), (5.20) 

where h > and t are fixed. From the Taylor theorem with integral remain- 
der, we have 

rt+h 



A h f(t) = [ f"(r)m h (T,t)dr, (5.21) 
Jt-h 
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where m/ l (r, i) is the hat-shaped function 

< + \ ( h-\r-t\, \r-t\<h , , 

m *( T ' t ) = \0, otherwise. ^ 

It is clear that mh(r, t) > for r £ [t — /i, t + /i]. Now, for any > 0, s > 0, 
1 > /i > 0, and t > t c + 1 + e, 

rt+h 

A h F N (t)= / F^(r)m h (T,t)dT 
Jt-h 
rt+h 

^ N (T)m h {T,t)dT (5.23) 



t-h 

K , t +h 



U-nJt-h iv 



2K-2\ 



k=0 ' 

where the error term is uniform in t > t c + 1 + e, since the asymptotic 
expansion f|5.16j> is uniform in t. On the other hand, 

A ^(t) = £ + 0(N^<-% (5.24) 

fc=0 

where the error term is uniform in t > t c + 1 + £, because the asymptotic 
expansion (|5.8|) is uniform in i. Since the coefficients of an asymptotic series 
in powers of N~ 2 for the function A^F^^t) are uniquely determined, we 
obtain from (j5T23l) . (ET241 that for any k > 0, 

£l2k(T)mh(T,t)dT. (5.25) 



t+ft i2fp{2k)( \ 

A h F^ k \t) = \H mh ( T ,t)dr, (5.26) 



But by (I5T2TT) . 

hence 

rt+h £,p(2k)u.\ ft+h 



dr 2 



m h {T,t)dr = g 2 k(r)m h (T,t)dT, (5.27) 



t-/i Jt-h 

and since h and i are arbitrary and m^r, i) is a positive function, equation 
(j5TT9]) holds. 

By Corollary 14.41 the function g 2 k( u ) is an even function of u, analytic at 
u = 0, hence as it — » 0, 

52fc («) = oo(*0 + a 2 (A:)n 2 + 0(u 4 ). (5.28) 

Because of (|5.16p . we know that 

= ^>0, (5.29) 
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and this implies that as t — > oo, 

+ W + ^ (5 - 30) 

Integrating twice we obtain that 

pm(t) = M*)^ _ aj^a + c + m + R2t((h (5 3l) 

where 

#2fe(t) = 0(t~ 3/2 ) (5.32) 
and C, -D are some unknown constants. The error term can be written as 

^w=££(»w-f$ («■») 

Combining ()5.31j) for = with ()5.9[) we obtain that 
F(o)(t) = ^-^M + F(o)(„) 



3 4 
4a (0)t 3 / 2 ai(0)lnt 



(5.34) 



+ C + Dt + R (t). 



3 72 
Since 

lim F^fu) = 0, lim R (t) = 0, (5.35) 

u — >0 t— >oo 

we conclude from (|5.34p . by taking t — > oo, that 

ao(0) = ^, oi(0) = 18, C = -^, D = 0, (5.36) 

and (|5TT7) follows from (|fT33jh (jOIj) . 

Similarly, combining (|5.3ip for > 1 with (|5.10p we obtain that 

F' 2k \t) = F< 2 ')( U ) 

= M^!_M^ +C+D(+B2i(t) . (5 - 37) 

Again, since 

lim F {2k) (u) = 0, lim R 2k (t) = 0, (5.38) 

u— >0 t— >oo 

we conclude from (|5.37p . by taking t — >• oo, that 

a (/c)=0, oi(jfe) = 0, C = 0, D = 0, (5.39) 

and ([EISD follows from (jQ3|) . ([5737]) . □ 

Remark 5.2. From (|5.28p and (|5,39p we obtain that if k > 1 i/ien 

52fc (n) = 0(n 4 ), u^O. (5.40) 

Also, from (|5.9p . (|5.10p . (|5.17p . and (|5.18p . we have the following corollary 
of Proposition 15.11 



28 



PAVEL BLEHER AND ALFREDO DEANO 



Corollary 5.2. 

^<">-££("'>-^-5*)**- (541) 

and 

Using this result, we will prove Theorem 11,21 by integrating the explicit 
expression that we found for the leading coefficient go(u). Also, we will 
prove Theorem 1 1.31 by using an expression for g2(w) obtained from the string 
equations. 

6. Proof of Theorem 11.21 and Theorem 11.31 
Observe that from (|5.3I) we have that 

72u 2 = t~ 3/2 , (6.1) 

so (|4.26p implies that 

2T(j + l)r(i^) ~p o 2T{j 
Multiplying by the factor (3u) 2 / 3 = we obtain 

i i 00 r(^±i)t^^ 
* w = i7J + 4? + giF(7tT)?(ffl)- (6 ' 3) 

Integrating twice in t, according to formula ()5.41f) . we get 

F(°)(u) = V 1( 2 jt _ (6.4) 

^(3i-i)(3j-3)r(j + i)r(^)' 1 1 

and making the change of variable t = 47^173 > we have 

. 3j-3 



^ 2 r(i + i)r(^) ^ o 2 r(j + i)r(^) 



00 3 2(i-l) 4 ^ rf 3i+lN 2(j-l) 

F (o)( u ) = y 4 2 n 2 ju 

^(3j-i)(3j-3)r(j + i)r(^ ; 

"^2r(j + 3)r(i + i)' 

This proves formula (jl.18)) . 

L'h ill ill I 1 "1 / 1 /""-li / 1 f I 1 / ■ 1 ,'\ 1 "1 i t.' ^ 



(6.5) 



From the explicit expression of the coefficients in formula (|1.18p we 
obtain the asymptotic formula, 



TOPOLOGICAL EXPANSION IN THE CUBIC RANDOM MATRIX MODEL 29 

which gives (I1.19D . together with some higher order corrections. Theorem 
11.21 is proved. 

When k = 1, system (|4.44p - (j4.45|) reads 

6 66oH - l\ (g 2 (w)\ = 3 / -5o'H \ fg7 x 

1-eSoH -630H y V^W A&M&gH/ ' l ' J 

This system can be solved (and simplified, using the string equations) to 
give 

. f \ 162goH(5 - 324g ( W )) 

52 M = 4 , (6.8 

(1 - IO850 {w)Y 

i 54w 

h ^ w) = gowa-mihwr (6 - 9) 

(6.10) 

An explicit expression f|l .21 [) for the coefficients fyj can be now obtained 
from equation (|6.8p . Namely, g2{w) can be expanded in powers of w around 
the origin, 

oo 

faiw) = J2 c j wJ i ( 6 - n ) 
i=2 

starting with a tu 2 term. From cubic equation (|4.39|) . we obtain 



w = g (w)y/l-72g (w), (6.12) 
so writing v = go(w) and applying the Cauchy integral formula, we find 

1 / (5-324u)«-J'(l - 72v) =s ^ L ( . — 36u \ , 

' ' Vl - 72w dv, 



2vri/ 7 162(1 - 108u) 4 V ~ VT^72v 

1 / (5 - 324w)«-J'(l - 72u)"2 



_ i 



2vrf J 7 162(1 - 108v) 3 

(6.13) 

where 7 is a smooth closed contour around the origin in the v plane. Now 
we expand the binomial series and pick up the residue at v = 0, to get 

j' — l /3j -I \ / O \ m 

Cj = 162 • 72^ ^ ( T J m ) ( m + i)( m + 5) ( £ ) . (6.14) 
m=0 V J - m - 1 / V - / 



If we write this in terms of the standard Pochhammer symbol, see pQ, we 
see that 

(rn+l)(m + 5)=5^4§ !!1 , (6.15) 



and also 

, , , ■ , , / ( 6 - 16 ) 
j-m-l \j - I J (_3i + ^ 
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so we can identify Cj with the following hypergeometric function: 
45.72^f-l^ 

' 3-^2 



j-1 



3 + 1,2,6 3 
22 + 1 '2 



5, -j 



.17) 



Since one of the parameters in the numerator equals another one in the 
denominator plus one, we can simplify this 3F2 function in terms of Gauss 
hypergeometric functions using the following straightforward identity: 



. a,b,c+l 
3i?2 ( c,d > Z 



2F1 



a,b \ abz 



a + 1,6 + 1 
d + 1 



and thus 



45_ ,/¥ - 1 



-72 J 



Vj-i 

j + 1,2 3\ , 6(j-l) 



-22 + 1 
2 ^ L 



+ 



2 1 ' 5(3j-2) 
Next, observe from (|4,35p and f|4.36j) that 



1F1 



-j + 2,3 3 
-| + 2 '2 



.18) 



.19) 



1 



2\/i 



52 («), 52 (u) = u 2 g 2 (w), 



.20) 



so 



h(t) 



1 x 

144t2 S C J 



1 



72t 3 /2 



72t 3 / 2 



(6.21) 



Using formula (|5.42p . we integrate twice in t, and going back to the variable 
u, we find 

3- 

so the coefficients are 



36 3j(3i + 2 ) 



(6.22) 



f (2) 
'2j 



ci(2j)! 



108j(3j + 2) 
45 • 72^ _ r 
J-1 
5-72?Tf¥)(2i) 



48(3j + 2)r(j + l)r 2 + 1 



-j + 1,2,6 3\ (2j)! 
5,-f + l '2/ 108j(3j + 2) 



(6.23) 



3-^2 



-J + 1,2,6 3 
5,-| + l '2r 



This proves formula (jl.2ip . The large j asymptotic behavior of these coef- 
ficients is given in Appendix ICl below, proving formula (jl.22p . 
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7. Case of genus g > 1 

The expressions for the terms in the topological expansion become more 
and more complicated as k grows. Fortunately, as we pointed out before, 
in order to understand the large j behavior of the terms f^ 9 \ it is enough 
to consider the leading order terms near the critical point w = w c only. We 
recall that w c = ^f, as given in f[4.40j) . 

When w = w c , the cubic equation 72gjj(u;) — g~o(w) + w 2 = given in 
(|4.39p has three roots, 

ri = ro = , ra = . (7.1) 

108 216 y ' 

Substituting 

w = w c -Aw, g (w) = + Ag (w) (7.2) 
into the cubic equation, we have 

72(Ag (w)) 2 (goW + 2jg) = w c ~ w2 = + w), ( 7 - 3 ) 

Therefore, 



(7.4) 



/ (w + w c )Aw / 2w c Aw m , . , 

Ago = ±i \ , r — i-T = ±W — / , \ rr + O Aw 

ol/2ol/4 

= ± (Aw) 1 / 2 + O(Aw) 

18 

So, near the critical value w = w c , our solution satisfies 

9oH = -L - ^(Aw) 1 / 2 + O(Aw), (7.5) 

taking the minus sign since the function go(w) is increasing with w and 
therefore go(w) < g(w c ). 

Since we have bo(w) explicitly in terms of go (id), see formula f|4.38j) . we 
can obtain the asymptotic behavior of bo(w) as well: 

So(«j) = ^— ^ ~ 2 1 / 2 3" 1 / 4 (Aw) 1 / 2 + O(Aw). (7.6) 
18 

We also note that near the critical point w = w c , the determinant D(w) in 
(I4.47P behaves as follows: 

D(w) = 2 3 / 2 3 5/4 (Aw) 1 / 2 + 0{Aw) = D'{Aw) 1/2 + 0{Aw). (7.7) 

We can write formulae (|7.5p . (|7.6p as 

= soK) + Cq(Aw) 1 / 2 + O(Ato), 

SoH = 6oK) + Do(Aw) 1 / 2 + O(Aw), 
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where 



9o(w c ) = — , b (w c ) = , (7.9) 



and 

18 



ol/2ql/4 

Co = , = -2 1 / 2 3- 1 / 4 . (7.10) 



Observe that 

D = 6V3C . (7.11) 
For higher order terms, k > 1, we make the following Ansatz: 

92k{w) = C 2k (Aw)^-^ + OdAw) 1 '^ ), 



(7.12) 



where the functions g 2 k(w), b 2 k{w) have a square root singularity in the com- 
plex plane at the point w = w c . We will prove this Ansatz by induction with 
respect to k, by using the system of equations (|4.44p - (|4.45p . Simultaneously, 
we will derive recurrence equations for C 2 k, D 2 k- 

Let us analyze the first sum on the right hand side of (j4.44jl . By differ- 
entiating the first formula in (I7.12p . we obtain 

= cSi\Aw) l 2~^ + ..., m < k — 1, (7.13) 
where the dots indicate higher order terms with respect to Aw and 

cffi> = ft.(^-i)...(^ + W -l)-i). (T.U, 
Since m + j = k, we can write formula (|7.13p as 

<?£?H = ^(A^l-f +1 + . . . , j > 1. (7.15) 

Hence the leading term in the first sum on the right hand side of (|4.44|) 
corresponds to j = 1, m = k — 1, so 

_ fi V 02m M 

^ (2j)!2 2 -? 

m+j=k; m<fc— 1 

= _^(^_l)(^ + l) (Aw)1 -f + ... (7.16) 

3C 2fc _ 2 (5fc - 6) (5fc - 4) u-f , 

16 ( 2 ... 

In the second sum, we have 

5fc 



^ b 2m (w)b 2m > (w) = -3 ^ D 2m D 2m i(Aw) 1 "2+... 

(7.17) 



m+m'=fe m+m'=fc 
m,m'<fc— 1 m,m'<k— 1 
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For k = 1 this sum is equal to 0. Substituting formulae (I7.16D . (I7.17D into 
(j4.44|) . we obtain 



<*■ a i 5fc 

(u>) - (1 - 66 (w))M™) = ^(Au;) 1 -- + . . . , 



(7.18) 



where 



. 3C 2fc _ 2 (5fc - 6) (5* - 4) 

^2fc = 77 -3 

16 



D 2m D 2m ,. (7.19) 

m+m'=fc; m,m'<k— 1 



Similarly, in equation (I4.45D the leading singular term on the right, with 

(Aw) 1_ ~, corresponds to j = 1, m = 0, m' = A; — 1 and to j = 0, m + m' = 
A;, so that 



(1 - 6b (w))g 2k (w) - 6g (w)b 2k (w) = B^Aw) 1 2 



(7.20) 



where 



£ 2fc _ 2 (5fc - 6) (5fc - 4) 
* 2fc = 576 + 6 



c 2m D 2m i. (7.21) 

m+m'=fc; m,m' <k— 1 



Solving system (I7.18|) . (|7.20j) and using asymptotic formula (I7.7j) for its 
determinant, we obtain that 



92k(w) = C 2k (Aw)2-- + ^((Au;) 1 --) 
hkH = D 2k (Aw)^ + Oi(Aw) 1 -^ ) 



(7.22) 



where 



C 2 k 
D 2k 



1 



2 3/2 3 5/4 
1 

2 3 /2 3 5/4 



6<7o(w c )^2& + (1 - 66o(^ c )) J B 2fc 
(l-6SoK))A 2 fc + 6B2fc 



(7.23) 



This proves Ansatz (|7.12p and gives the recurrence formula for the coeffi- 
cients C 2k , D 2k . Moreover, substituting formulae (|7.9p . we obtain that 



C 2k 
D 2k 



1 



A 2k V3B 2k \ 
2 :s V " 18 3 / 

I / V3A 2k 



2 3 /235/4 



+ 6B 



2k 



(7.24) 



Hence 



D 2fc = 6V3C 2k , 



(7.25) 
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which extends relation (17. lip to A; > 1. By using this relation, we can 
eliminate D 2m from formulae (]7.19p . f|7.21[) : 

3C 2fc _ 2 (5fc-6) (5fc-4) ^ 

A 2k — 2-~i ^2m^2m', 

m+m' =k; m,m'<k—l 

V3C 2k „ 2 (5k-6)(5k-4) Qc f- 

£>2k — gg hOOVO L'2m^2m', 

m+m'=fc; m,m'<fc— 1 

(7.26) 

and hence by the first formula in f|7.24j) . 

C ^ = ¥ J 2 ^((5k-6)(5k-4)^^ + 54 Y, C 2m C 2m \ 

m+m'=k 
m,m' <k— 1 

(7.27) 

This proves recursive formula (|1.27p . 

From (|7.10p and (|1.27p . we find the values of the constant C 2k for k = 
0,1,2 as 

2V^ 1 49 • 2V^4 
18 - ' C2 "48Tl08' Q " 17915904 " 1 j 

Let us apply formula (j7.12j) to find the asymptotic behavior of F^ 2k \u) 
at u = u c for k > 1. From (14.36P and (I5.29p . we know that 

~ M 92k{u) u 4k ~ 2 g 2k (w) mtg 2k {w) 

92k{t) = ^Vf= 2Vt = (7 - 29) 

and then we obtain from (17.221) that 

hk(t) = C 2k (Aty2-f + ..., (7.30) 

where 



a 



2k 



1 5k 

36t c C 2k ( dw \2~ — 



t=t. 



{72Hl) k V dt 



(7.31) 



If we integrate (|7.30p twice with respect to At, we get 



if k > 2. By (I5.10p we obtain now that 

F m {u) = A2k (^_^y 2 +..., k >2, (7.33) 



where 



4C 2k u 5 c ~ 5k ( dt \5"f fc 

w=w c 



A - - 



(5Jfe-5)(5fc-3) V <i«> 



(7.34) 
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Observe that by (16. ip . 

1 _ 2 /3 dw 

W= 72t^ 1 t c = 3-2 , - — 



t=t n 



48t 



5/2 



hence 

144 • 48 2 tS^ / 1 x k 



u cJ ^ u; 



2j 

3=0 



where c,- has the following asymptotic behavior as j — > oo: 



c 



j(5k-7)/2 
i ~ y f5k-5\ 



2 / 

Combining (03]) . (|7^71) ~ (I7\39| ). we find 



This proves Theorem 11.51 
The initial values of K 2 k are 



(7.35) 



A2k= (5k-5)(5k-3) {WMJ ^ (7 ' 36) 
Substituting expression (|7.3ip and simplifying, we obtain 

_ 24 • 3^C 2k 

Mk ~ (5k-5)(5k-3)u* ' {7 - 67) 

Using the binomial expansion, we have 

7 ,2\f-| fc °° c . u 2j 



;i+C(r 1 )). (7.39) 



oo 2 7 

F^)(n) = j:4f\, (7-40) 
• 1 w c 

3 = 1 L 



4f = W 5fe - 7)/2 (i+o(r 1/2 )), (7.4i) 



where as j — ^ oo, 
with 

6 • 3 1 / 4 C 2fc , v 



= -i= , #2 = 4r , K± = ^-== . (7.43) 

V6t7 48' 1440^/oV V ' 



When k = 1, formula (|7.33p becomes 
where 

C3 = ^-(l + 0(r 1/2 )) • (7.45) 
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8. Counting 3-valent graphs on a Riemannian surface 

In this final section, we prove that the coefficient /i in series (|1.17p is 
equal to the number of 3-valent graphs with p vertices on a closed Riemann- 
ian surface of genus g. By differentiating formula (|1.5|) p times with respect 
to u and evaluating at the point u = 0, we obtain 

where 



vi=l 



{f(z 1 ,...,z N )) = z I ... I f{zi,...,z N ) 



N ^ (8.2) 

x n fa ~ Zfe ) 2 n e 2 dzi ■ ■ ■ dzN 

is the mathematical expectation with respect to the Gaussian ensemble. By 
the Cauchy theorem, we can deform the contour of integration V in the 
integral on the right in (|8.1|) to the real axis, and then we can return to the 
matrix integral: 



/ iyi\ \ =((TrM 3 f) =^i- / (TtM 3 ye- NT ^dM. 

(8.3) 

For odd p the latter integral is equal to 0. To evaluate the integral for even 
p = 2q, we apply the Wick theorem. 
We have 

N 

TrM 3 = M i3 M jk M ki (8.4) 

i,j,k=l 

and 

N 

(TrM 3 f= J2 M iljl M jlkl M klil ...M ipjp M jpkp M kpip . (8.5) 

i-lljl 5 ■ ■ ■ }Jp )^p — 1 

The covariance matrix of Mjj's is 

(Mi3M kl ) = ^* . (8.6) 

By the Wick theorem, the expectation of the product of matrix entries with 
respect to the Gaussian measure can be expressed in terms of the product 
of expectations of all possible pairings of the matrix entries (see for instance 
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[2lT §1.6] for a more detailed exposition). In this case, 
(M hjl M jlkl M klh . . . M ipjp M jpkp M kpip ) Q 

3p/2 3p/2 /g rj-j 

= £ II (M iis)j ( s) M k(s)lis) ) o = Y^N-^' 2 H 5 iis)lis) S jis)kis) . 

7T S = 1 7T S=l 

where the sum is taken over all partitions 

vr : {M hjl , M jlkl , M klh M ipjp ,M jpkp , M kpip } 

^ (8.8) 
= \J{ M ^j^> M k^i(°)} 

s=l 

of the set {M hjl , M jlkl , M klh , . . . , M ipjp , M jpkp , M kpip } into disjoint pairs 
{ M H")j(s), M k{s)l{s) }. 

From (18751) . {I57TD we have that 



N 3p/2 
n h,jiM>—,ipdp,k p =l s=l (8.9) 

7T 

where / is the number of cycles in the set 

{Midi > M nk, > , . . . , M ipip , M, pfcp , M fcpip } 
generated by the equalities 

i(')=l(') t j( s )=k^ s \ 8 = l,...,3p/2. (8.10) 

Thus, 



Z N (0) ^ ' 2 

By the second Wick theorem, we obtain now that 

F - )(o)= ^ in €^ = ? v+/ "" (s - i2) 

where the sum is taken over partitions tt such that the graph r(7r) is con- 
nected (for a general formulation of this result in combinatorics see for in- 
stance [Ml Vol. 2]). 

From (TLT6j) . ([07]) we have 

*o)-EC' (8 - 13) 

Comparing this with (|8,12p . we obtain that for each partition ir, 

p + f-l-2 = -2g (8.14) 
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Figure 3. 3-valent connected ribbon graphs with two ver- 
tices, genus 0. 



and 

4 29) = E c !> ( 8 - 15 ) 

tt: p+f-l-2=-2g 

As was noticed in [7J, by the Euler formula, the number g in (|8.14j) is equal 
to the minimal genus of a closed oriented Riemannian manifold on which 
the graph r(7r) can be realized without self- intersections. See the work |23| 

of Mulase for a rigorous proof. This proves that fp 2g ^ counts the number of 
3-valent graphs with p vertices on a closed oriented Riemannian manifold of 
genus g. 

As an illustration, Fig. 2 depict possible 3-valent (ribbon) graphs with two 
vertices on the sphere, and Fig. 3, on the torus. Observe that by counting 
the possible enumerations of the edges, we obtain the multiplicity factor of 
the first graph in Fig.2 to be equal to 9, and of the second graph to be 
equal to 3. The total number of the graphs with multiplicities is equal to 
12, which fits well to the first coefficient 6 = in the expansion of f(0)( u) 
given in (|1.20p . Similarly, we obtain that the multiplicity factor of the graph 
in Fig. 3 on the torus is equal to 3, which fits well to the first coefficient | 
in the expansion of F( 2 \u) presented in (fQ3|) . 
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Figure 4. A 3-valent connected fat graph with two vertices, 
genus 1. 



Appendix A. RHP and properties of the orthogonal 

POLYNOMIALS P n {z) 

In this appendix we prove various propositions from Section [3j 
Proof of Proposition \3.1\ Applying equations (j3.3H to the last column in 
the determinant for D n , we obtain that 



co 
ci 

C n -1 

co 
ci 

Cn-1 
C/t 



ci 

C2 

( '-n 
C n +1 

Cl 
C2 

Ctj 

c n +i 



Cn-1 
Cn 

C2n-2 
C2n-1 

Cn-1 
C?f 



C2n-2 
C2n-1 



C/i 
Cn+1 

C2n-1 
C2n 






/?.„ 



(A.l) 



n—li 



where 



which proves (13.8 



z n P n (z)e- NV ^dz = / P n (z) 2 e - W ^cte, 



(A.2) 



□ 



Proof of Proposition 1 3. & Let y n (z) = det Then by (|2.28jl . y n+ (z) = 

y„_(z), z £ Tq UTi, hence y n (z) is an entire function. By (|2.29j) . y n (z) — > 1 
as z — > oo. Hence, by the Liouville theorem, y n {z) = 1. 

Next, we claim that Y n (z) is the unique solution of the RHP. Indeed, 
suppose that Y n ' 2 (z) are two solutions. Then X{z) = Y^ l {z)Y^{z)~ 1 has 
no jump on T\ U T2 and X(z) — > I as z — > oo. By the Liouville theorem, 
X(z) = 1, hence Y^(z) = Y^(z), so the solution is unique. 
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Consider now the (ll)-element Y n n{z) of the matrix Y n {z). From (|2.28p 
we obtain that Y n u+(z) = Y n u-(z), hence Y n \\(z) is an entire function. 
By fZlOD , Y nU (z) = z n + O^ 1 - 1 ), hence Y nn (z) is a monic polynomial of 
degree n. 

Consider then the (12)-element Y n \ 2 {z) of the matrix Y n (z). From (|2,28|) 
we obtain that for z £ Tq U Ti, 

Y nl2+ (z) - Y nl2 ^(z) = a(z)e- NV ^Y nll (z), (A.3) 

and from (|2.30p . that as z — > oo, 

Y nl2 (z) = 0{z- n ~ l ). (A.4) 

This implies that 

i r e _w ( s )y 11 (V) 

2m J T s — z 

where again r = To U Y\. By expanding l/(s — z) into a geometric series, 
we obtain that 

Comparing this with (|A.4p . we obtain that 

J s k e- NV( - s) Y nll (s)ds = 0, k = 0,l,...,n-l, (A.7) 

hence Y n n(z) is a monic orthogonal polynomial of degree n, Y n u(z) = P n (z). 
Since the solution Y n (z) to the RHP (I2.28|) - (l2.30p is unique, the orthogonal 
polynomial P n (z) is unique as well. □ 
Proof of Proposition \3.3[ Consider the matrix- valued function 

U n {z) = Y n+l (z)Y n (z)-\ (A.8) 

It has no jump, hence it is an entire function. As z — > oo, U n (z) = O(z), 
hence U n (z) is a linear function. Moreover, (]2.30p implies that U n (z) has 
the form 

v M -(' + ™ (A,) 

Considering the equation Y n+ \(z) = U n (z)Y n (z) for the (l,l)-entry, we ob- 
tain the three term recurrence relation (|3.1ip with some coefficients, which 
we denote j3 n and r y 2 . □ 
Proof of Proposition \3.4\ Integrating by parts, we obtain that 



= / P n (z)P' n (z)e- NV(z) dz = N P n {z) 2 V'(z)e- NV(z) dz 

Jr Jr (A.10) 

= N P n (z) 2 (z - 3uz 2 )e- NV{z) dz. 
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Now, from (13. lip we have that 

J P n {zfze- NV ^dz = h n p n , (A.11) 

and 

P n (zfz 2 e- NV ^dz = I [P n+l (z)+p n P n {z) +1 lP n ^{z)} 2 e- NV ^dz 
Jr 

= h n+ i + 0ih n + ^nK, 

hence 



(A.12) 



= h n (3 n - 3u(h n+1 + fihn + 7X-i)- (A.13) 

This implies the first equation in (|3.14p . 

Similarly, integrating by parts, we obtain that 

nh n ^ = J^P n ^{z)P' n (z)e- NV ^dz 

= N J P n ^ 1 (z)P n (z)V'(z)e- NV ^dz (A.14) 

= N j P n _!(z)P n (z)(z - Zuz 2 )e- NV ^dz. 

Now, from (|3.1ip we have that 

J P n ^(z)P n (z)ze' NV{z) dz = h n , (A.15) 

and 

j^P n ^(z)P n (z)z 2 e- NV ^dz 

= j[P n { Z ) + Pn-lPn-l(z) +lLlPn-2(z)]Pn(z)ze- NV ^dz (A - 16) 
= {Pn + Pn-l)h m 
hence 

nh n -i = N[h n - MPn + Pn-i)h n ]. (A.17) 

This implies the second equation in ([3. 14)1 ■ □ 
Proof of Proposition \3. 5[ Observe that the moments, 

z j e -N(* 2 /2-™ 3 ) dZj (A.18) 

are C°°-functions of u for u > 0, and they are analytic for u > 0. Hence, 
the same is true for the determinant P> n -\ and the polynomial D n (z), for 
any n. □ 
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Appendix B. Proof of Toda equation (15. ip 
Let N be fixed. Then, as u — > 0, the complex- valued measure 

e- N ^~ uzZ Uz 

on r = aTo + (1 — a)ri converges to the Gaussian measure, and the mo- 
ments of the measure converge to the moments of the Gaussian measure. 
This implies that there exists u(N) > such that for u G [0, u(N)] the or- 
thogonal polynomials P n (z) exist for n = 0, 1, . . . ,N. Hence the orthogonal 
polynomials P n (C) exist for n = 0, 1, . . . , N, if t > t(N), where 

t(N) -li^W (R1) 

cf. (|5.10p This implies, by the usual argument (see e.g. {5]), that Toda 
equation (|5.15p is valid for t > t(N). Since both the free energy Fn and 
the recurrence coefficient 7^ are analytic in t for t > t c , we obtain, by the 
analytic continuation, that Toda equation (|5.15p is valid for t > t c . □ 



Appendix C. Proof of formula (|1.22p 

In order to prove (|1.22p , we will use formula ()6.18j) and some linear trans- 
formations and integral representations of the 2-^1 functions. Namely, we 
note that we can use [TJ 15.8.7]: 



2^1 



-m,b \ _ (c - b) m F ( -m, b _ 

C ' J (c) rn \b-C-m + l' 



setting m = j — 1 and m = j — 2: 



-j + 1,2 .3\ V^ + Vi / -j + 1,2 , 1 
¥ + 1 '2j-(I7^ 2i?1 U + 3 '"2.) (CJ) 



and 



-j + 2,3 S\ _ V2 + V2 / -J + 2,3 , 1 
2 l( -1 + 2 '2j"/ 1 + 1 \ 2jFl U + 4 ~2j- (( -' 

V 2 



Additionally, we can use 
1,6 \ 

c' Z J T(b)T(c-b)J {1-ztY 
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which is valid provided that Re c > Re b > 0. In our case, 

-J + 1,2 3\ r (l + 3) yi ^ , / l\H 



F i(i,o' ;o) = — h — f / t(i-i)5(i + ^l 

(C.4) 



x 2 

If we write the previous integrands as follows, 

h = + Itj 1 e-^dt, h = f Q t 2 (l + ±t\ ^ e-^dt, 

(C.5) 

with 

<f>(t) = -^ln(l-t) -In (l + lt J . (C.6) 



2 V 2 , 

This function has a minimum at £ = 0, since 

0'm = _ — ii — _ (c.7) 

w 2(1 — *)(2 + *)' V ' 

which is positive for < t < 1. Therefore, for large j the main contribution 
to both integrals comes from t = 0. We make the change of variable 

--ln(l-i)-ln(l + -H =t, (C.8) 



which maps t = to r = 0. We expand around t = 0: 

_ I ln( i _ _ hl (i + ^ = |ta + ^3 + ^4 + 0(t 5 }) (C 9) 

so inverting this series, we obtain 

t = ^ T V* --t- ^r 3 / 2 + 0(r 2 ). (CIO) 
3 9 81 v ; v ; 

Additionally, one can work out the following: 

r 3/2^ 



2 / t\ 2 dt 8\/6 , /a 160 
V 2/ dr 9 27 



(C.ll) 



T 



An application of classical Watson's lemma [23] , shows now that the integrals 
in (|C,5|) are I\ = 0{j~ l ) and I2 = 0(j~ 3 / 2 ) respectively. The factors that 
multiply in front are 
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and 

Ci + 4 X- 2 r(| + 4) _r(f + 2 ) _ 27j3 



i + i) 3 ._ a 2r ( i + 1) 2r (a _ i) 16 

so we only need to consider the second integral I2, which is dominant for 
large j. Writing everything together and noting that the binomial number 
has the following behavior, 



j- 1 ) v{j)vU + i) 



, 2 

we obtain 



j -> 00, (C.14) 



45-72' r (,Tj 6(j-l) „ / -i + 2,3 _3 



r(i)r 1 + 1 



when j — t- 00. We observe that we can write this as 



(C.15) 



2- 2 3 3 j(3^2 2 ) J ' = ^l, (C.16) 



and recalling (|6.22p we obtain 



g, ^ 27i(2j)! _ 
J 4 • 36 • 9j 2 ^ 48j^ 



which proves asymptotic formula (jl,22p . 
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